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We have calculated high temperature series to 12th order in inverse temperature for singlet su- 
perconducting correlation functions of the 2D t-J model with s-, d^2_y2- and dai^-symmetry pairs. 
We find the correlation length for d^2_y2 pairing grows strongly with decreasing temperature and 
develops a broad peak as a function of doping at T / J — 0.25 for J/t = 0.4. The correlation lengths 
for s- and d^^-symmetry remain small and do not display peaks. Antiferromagnetic spin correlations 
at low doping act to suppress the d^2_y^ and dxy superconducting correlation lengths. Our results 
support the hypothesis that the strong electronic correlations found in the Cu02 planes of high 
temperature superconductors are the origin of the superconducting order. 



Shortly after the discovery of high temperature super- 
conductors by Bednorz and MiillerQ, Anderson pro- 
posed that the new superconductors should be viewed as 
single band doped Mott insulators, with a novel, non- 
phonon origin for their superconductivity. The model 
used to describe the Cu02 planes found in these ma- 
terials is the two-dimensional Hubbard model, which in 
its strongly correlated limit was restated by Zhang and 
RiceQ as the 2D t-J model. Below we show that the 2D 
t-J model has strongly enhanced d,j.2_y2 pairing correla- 
tions for physical choices of the model parameters and 
doping range. 

The Hubbard and t-J models have been widely studied 
by every method available in theorists' arsenals^. The 
results to date for pairing co rrelations have been at best 
mixed 0, E S i, ES Enil m ■ Lacking small parame- 
ters, strongly correlated models present severe difficulties 
for analytic calculations and numerical techniques have 
difficulties with finite size effects and the fermion sign 
problem. To avoid these problems, we calculate high tem- 
perature series to 12th order in inverse temperature /3 for 
the pairing correlations in the 2D t-J model. The series 
are generated by a linked cluster expansion [l^. putting 
four extra operators on each cluster for the two pairs. 

The Hamiltonian for the t-J model is 



and the Hermitian adjoint 



H = -tP 



(ij>o" 



(4c,.+4c,,)p + J^S,.S„ (1) 



where the sums are over nearest neighbor pairs of sites 
on a square lattice and the P operators eliminate states 
with doubly occupied sites from the Hilbert space. The 
spin operators are given by = S^^j cI^cTq^Cj^, where 
(Tap is a vector of Pauli matrices. 

We limit our calculation to spin singlet pairing, with 
the pair operator for sites i and j given by 



n/2 



(2) 



AT. 



(3) 



The Hermitian real space four-point correlation func- 
tion with operators on sites i, j, k and I is 



P{id,k,l) = i(Al^.A,i-f A|,A,,) (4) 

(4iCU>(c]|Cfet) -t- (4iC/ci)(c]|Q|) 
+ {4]Cll){c]^Ckl) + (4|Cfct)(c]^Qx) 

+ (c/\c,x)(4^Cj|) + (4|Ci|)(c|^Cjx) 

+ (c|tC»T>(4iCji) + (4iCii>(c|tCiT)J ' 

where we have removed the single particle contributions 
to P(i, j, fc, I) to focus on pairing fluctuations. For a non- 
interacting system P{i,j, k, I) — 0. Due to the no double 
occupancy constraint in the t-J model we have two re- 
strictions on the site choices: i j and k ^ I. Otherwise 
the pairs are allowed to take on all arrangements and 
sizes on the clusters. 

From our data we calculate series for the real space 
correlators P(r, r',R), where r and r' are vectors giv- 
ing the orientation and internal size of the two pairs and 
R is the distance between the pairs, determined by the 
distance between the centers of mass. Since the pairs 
are spin singlets, r and r' are limited to the upper half 
plane and positive x-axis to avoid double counting. All 
of the distances are measured in terms of the lattice spac- 
ing a. Since the original four point correlation function 
P{i,j, k, I) is identically zero for non-interacting systems, 
we also have P(r, r', R) = for non-interacting systems. 
The internal symmetry of the pairs is included by multi- 
plying P(r, r', R) by phase factors <j>i{r) and (piir') where 
each phase factor can have the values -fl, —1 or 0. We 
have (f>i{r) — when r is aligned with a node. Here the 
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subscript i labels the different symmetries we consider: s, 
dx2-y2 or dxy The phase factor of a single pair is shown 
in Fig. 1 for the anisotropic symmetries d^^-y^ and d^y- 
An s-symmetry pair has uniform phase, 0s (r) = +1 for 
all r. 

From P(r, r', R) and the phase factors we want to con- 
struct a correlator Pi(R) which depends only on R, the 
separation between the pairs, and the symmetry of the 
pairs {i — s, d.j;2_y2 or dxy)- This is then used to measure 
the strength of the pairing correlations by calculating the 
q = correlation length via 



1 ErIRP^.(R) 

2d Er^.(R) 



(5) 



with dimension d — 2. 

One way to calculate Pi(R) is to sum over r and r' 
independently 



P,(R) = 5]0,(r)^0,(r')P(r,r',R) 



(6) 



We find that for the anisotropic pair states with |R| > 
Eq. 6 gives Pi(R)'s that are very small and negative. 
This suggests there is spurious cancellation occurring in 
our definition of Pi(R). 

As an example of the problem that can arise when 
calculating Pi(R) by Eq. 6, consider the special case0 
of pairs limited in size to nearest neighbors with d^2_j,2 
symmetry. This is the type of calculation done in Refs. 
lilSllBlillSIIlIIland Ref. US We have two types 
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FIG. 1: Internal nodes and spatial phases for d-symmetry 
pairs: a) d^2_y2, b) d^y The light lines are the underlying- 
square lattice and the solid black lines are the nodal lines. To 
determine the phase, consider one operator of a pair at the 
origin (where the nodes cross) and put the second operator on 
a diflerent lattice site. The phase associated with that part of 
the lattice is the internal phase of the pair. The second oper- 
ator cannot be put on the nodal lines; the pair wave function 
is zero for those configurations. For s-symmetry there are 
no nodes determined by symmetry and the phase is always 
positive. 

of correlators to consider (omitting the subtracted single 
particle terms for clarity): 

(0+) = (c|^c|^xTC''+R-+XTC»+Ri) 

(0_) = (c|^c|^-^c,+R+j^|C,+Rl), 



where x and y are unit vectors in the x and y directions. 
For dx2_y2 symmetry we have the phases 0+ = -1-1 and 
= —1. Now consider what happens to Pa;2_j^2(R) 
above and below Tr as IRI — > oo: 



and 



> |$o|' + A<cxp(-|R|/e,2_y2) r<Tc 

> A>exp(-|R|/e,2_j,2) T>Tc 

-|$oP + 5<exp(-|R|/e,2_^2) r<re 

B>exp(-|R|/42_^2) T>Tc, 



where |$oP is the magnitude squared of the dj.2_j,2 su- 
perconducting condensate and the A's and P's are non- 
universal coefficients of the parts of the correlator decay- 
ing with correlation length S,x^_y2. Below Tc we then 
have P^2_y2(R) ~ (0+) - (0_) - |$oP as |R| ~* oo due 
to the compensating phases for (0-). This is the defi- 
nition of long range superconducting order with dx2-y2 
symmetry. 

High temperature series, however, are limited to the 
non-ordered high temperature phase where Px2 _y2 (TV) ~ 
(A> - P>) exp(-|R|/^3;2_j^2). Now if Ay and P> do 
not refiect the symmetry of the ordered state we should 
expect Px2_y2{'R) to be typically small and of indetermi- 
nate sign. Extracting S.x2-y2 via Eq. 5 is essentially im- 
possible. In our calculation, (0+) and (0-) are both pos- 
itive, giving exactly this problem. This difficulty should 
be expected if there is anything in a calculation that pre- 
vents the system from being in the long range ordered 
state. Eq. 6 is generally a poor choice for a correla- 
tion function outside of the ordered phase since it can be 
anomalously small for reasons unrelated to the strength 
of the correlations. 

To reduce the possibility of spurious cancellations in 
the definition of Pi(R) we choose to limit the P(r, r', R) 
we include to having r = r'. By doing this we hold the 
internal degrees of freedom of the pairs fixed while we 
measure the distance dependence R of the correlation 
function. We have one pairing operator with a definite 
internal configuration r at position zero and a second 
pairing operator with the same internal configuration at 
position R. We then sum over r since we are only in- 
terested in the R dependence. Calculating a correlator 
where we allow r and r' to vary independently mixes dis- 
tance correlation with internal degrees of freedom, which 
isn't the quantity we want to measure. Our definition of 
P«(R) is 



P.(R) = ^02(^)p(^^^^j^)^ 



(7) 



With this definition of Pi(R) we are guaranteed to not 
have cancellation in the r sum due to trivial symme- 
try reasons. The effects of different symmetries are still 
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FIG. 2: Examples of configurations of two pairs on a square 
lattice. The first pair is represented by x 's and the second by 
o's. a), b) and c) are examples of contributions to the R=(l, 
0) correlator. For s-symmetry a), b) and c) all contribute, 
while for d^2„y2-symmetry the nodes eliminate c) and only 
a) and b) contribute. For d^^y-symmetry the nodes eliminate 
a) and b) so only c) contributes, d) is an example of a possible 
pair correlation not included in our definition of Pi(R), Eq. 
7. 



present through the nodes of the pair wave function 
where (f>i{r) = 0. The correlator defined in Eq. 7 is 
more versatile than the correlator defined by Eq. 6 due 
to Eq. 7 being able to find growing fluctuations in the 
disordered phase and also show the order paramenter in 
the ordered phase. An important feature of Eq. 7 is that 
to distinguish s and (1^2 _y2 pairing we need to include 
pairs with |r| > 1. Also, to have S^'^y positive we need to 
include pairs with |r| > I. 

With PiCR.) defined by Eq. 7, series can be calculated 
for defined by Eq. 5. The resulting series are extrap- 
olated to low temperatures using Pade approximantsjl5| 
either in the original expansion variable (3J or in the 
transformed variable w = tanh(Q;/3J) used by Singh and 
Glenister(I6l]. Here a is a numerical factor which can be 
adjusted to improve the convergence of the Fades. The 
estimated errors are from the spread of the Fade approx- 
imants at low temperatures. 

Our results at T/ J = 0.25 are shown in Fig. 3. Clearly 
the 2D t-J model does not have s- or (ia;y-symmetry pair- 
ing for the range of parameters shown. However, the 
correlation length for (i2;2_j,2-symmetry pairing ^3.2_y2 is 
large and suggests the ground state of the 2D t-J model 
has long range superconducting order with (1^2 _y2 sym- 
metry. Further support for this interpretation is shown 
in Fig. 4, where £,x^-y2 versus doping is plotted for a 
range of temperatures. The temperature dependence 
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FIG. 3: Pairing correlation length versus doping. The t-J 
model parameter ratio is J jt = 0.4 and the temperature is 
T I J = 0.25. Results for s-, d^2_y2- and d^jz-symmetry pairing 
are shown. The lines are guides for the eye. 



shows ^x^-y2 grows rapidly with decreasing tempera- 
ture over most of the doping range shown in Fig. 4, 
though ^j.2_j,2 decreases with decreasing temperature for 
S ^ 0.04 This is distinct from s- and d^jy-symmetry pair- 
ing, which for low temperatures have and ^xy slowly 
decreasing with decreasing temperature over the whole 
doping range shown. The temperature scale for growth 
in ^3,2 _y2 also agrees well with the temperature scales 



found for two-point correlation functions 
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FIG. 4: Pairing correlation length for d2.2_y2-symmetry ver- 
sus doping for a range of temperatures. The t-J model pa- 
rameter ratio is J jt — 0.4. The data sets are in the order 
indicated on the plot, with the highest curve corresponding 
to T/J = 0.25. The lines are guides for the eye. 

Our results are in accord with experiments on high 
temperature superconductors^^. The pairing has 
dx2-y2 symmetry, with the doping range where S,x2-y2 is 
large comparable to where Tc is observed. The maximum 
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in £,x2-y2 is shifting to smaller doping as the tempera- 
ture is reduced, and likely is near i5 « 0.2 for T/J « 0.1, 
though calculations at lower temperatures are needed to 
show this. Also, ^a;2_j,2 grows very steeply on the under- 
doped side, while falling off more slowly for overdoping. 
This is due to ^ = at half filling for all symmetries 
and all temperatures and the strong antiferromagnetic 
correlations present for underdoping. In experiments Tc 
frequently jumps from zero to a non-zero value for un- 
derdoped samples with a small change in doping. 

Band structure estimates |21i] for t-J model parameters 
give a next nearest neighbor hopping term t' roughly the 
same size as J. This term has not been included in the 
calculations presented here. We expect that t' would 
allow holes to move more easily through the antiferro- 
magnetic background present in the t-J model, and thus 
would move the peak in S,x2-y2 to lower doping. Although 
a t' term is not necessary to have strong superconducting 
correlations, calculations including a t' term should im- 
prove detailed comparisons between experiments on high 
temperature superconductors and the t-J model. 

The 2D superconducting correlation length does not 
need to grow to an extreme size like £,af at half filling [l9| 
to produce true 3D long range order due to the effective 
interplanar hopping t± driving the 2D to 3D crossover 
for superconductors being much larger than the inter- 
planar spin coupling J±. The temperature dependence 
for ^2;2_j,2 shown in Fig. 4 is likely to produce Tc's in 
the range observed for high temperature superconduc- 
tors. Lower temperature results are needed to show this 
conclusively. 

A peak in S,x'2-y2 has not been observed in previous 
calculations. The most likely reason is the problem 
with cancellation in -PiXR) discussed above. Previous 
calculationsH SUHliliESEilEllli the authors 
have examined have mixed terms for anisotropic pair- 
ing correlators. The negative U Hubbard model does not 
have this problem due to the pairing being predominantly 
on-site and isotropic |23|. The q = pairing susceptibil- 
ity was calculated by high temperature series in Ref. 15, 
with the conclusion that dx2_y2 pairi ng i s not found in 
the 2D t-J model. The results of Ref. [13 are influenced 
by single particle terms we have removed and the cancel- 
lation effect. We believe these effects distort the results 
of Ref. and when properly accounted for the q = 
pairing susceptibility will show (1^2 ^y2 correlations are 
present in the 2D t-J model. 

In conclusion, we find strong (ij.2_y 2 -symmetry pairing 
correlations in the 2D t-J model for physical values of 
the model parameters and doping. This suggests the ori- 
gin of high temperature superconductivity is the strong 
electronic correlations found in the CUO2 planes. 
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Note added in proof. After submission of this paper we 
received Ref. |22 where calculations related to ours are 
reported. 
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